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Abstract 

In these lecture notes we present some connections between ran- 
dom matrices, the asymmetric exclusion process, random tilings. 
These three apparently unrelated objects have (sometimes) a similar 
mathematical structure, an interlacing structure, and the correlation 
functions are given in terms of a kernel. In the basic examples, the 
kernel is expressed in terms of orthogonal polynomials. 
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1 Gaussian Unitary Ensemble of random ma- 
trices (GUE) 

Initially studied by statisticians in the 20's-30's, random matrices are then 
introduced in nuclear physics in the 60's to describe the energy levels distri- 
bution of heavy nuclei. The reader interested in a short discussion on random 
matrices in physics can read [37], in mathematics |13], and a good reference 
for a more extended analysis is [38] . 

1.1 The Gaussian Ensembles of random matrices 

The Gaussian ensembles of random matrices have been introduced by physi- 
cists (Dyson, Wigner, ...) in the sixties to model statistical properties of the 
resonance spectrum of heavy nuclei. The energy levels of a quantum system 
are the eigenvalues of an Hamiltonian. For heavy nuclei some properties of 
the spectrum, like eigenvalues' spacing statistics, seemed to have some reg- 
ularity (e.g., repulsions) common for all heavy nuclei. In other words there 
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was some universal behavior. If such a universal behavior exists, then it has 
to be the same as the behavior of choosing a random Hamiltonian. Moreover, 
since the heavy nuclei have a lot of bound stated (i.e., where the eigenvalues 
with normalizable eigenfunctions), the idea is to approximate it by a large 
matrix with random entries. 

Even assuming the entries of the matrix random, to have a chance to 
describe the physical properties of heavy atoms, the matrices need to satisfy 
the intrinsic symmetries of the systems: 

1. a real symmetric matrix can, a priori, describe a system with time rever- 
sal symmetry and rotation invariance or integer magnetic momentum, 

2. a real quaternionic matrix (i.e., the basis are the Pauli matrices) can be 
used for time reversal symmetry and half-integer magnetic momentum, 

3. a complex hermitian matrix can describe a system which is not time 
reversal invariant (e.g., with external magnetic field). 

We now first focus on the complex hermitian matrices case, since it is the 
one we are going to discuss later on. 

Definition 1. The Gaussian Unitary Ensemble (GUE) of random ma- 
trices is a probability measure P on the set of N x N complex hermitian 
matrices, with density 

¥{dH) = ^ exp A Tr{H')^ dH, with (3 = 2, (1) 

where dH = YliLi ^Hi i Hk^jxat dRe(i7ij)dIm(ifj j) is the reference mea- 
sure, and Z]\f is the normalization constant. 

The meaning of /3 = 2 will be clear once we consider the induced measure 
on the eigenvalues. 

Concerning the meaning of GUE, notice that the measure ([1]) is invariant 
over the unitary transformations and has a Gaussian form. The invariance 
over the unitary transformations is physically motivated, because the physical 
system does not depend on the choice of basis used to describe it. Indeed, by 
imposing that the measure P is (a) invariant under the change of basis (in the 
present case, invariant under the action of the group of symmetry U{N)) and 
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a = 1/2N 


a = 1 


a = N 


Largest eigenvalue 
Eigenvalues density 


2N + C(iVi/3) 
Oil) 


C)(iVi/2) 


v^+C(Ar-2/3) 
0{N) 



Table 1: Typical scaling for the Gaussian Unitary Ensemble 



(b) the entries of the matrices as independent random variables (of course, 
up to the required symmetry), then one gets the family of measures 

exp (-aTr(/J^) + 6Tr(if) + c) , a>0,b,ceR. (2) 

The value of (c) is determined by the normalization requirement, while by 
an appropriate shift of the zeros of the energy (i.e., H H — E for some 
given E), we can set 6 = 0. The energy shift is irrelevant from the physics 
point of view, because by the first principle of thermodynamics the energy 
of a system is an extensive observable defined up to a constant. The value of 
a remains free to be chosen. Different choices can be easily compared since 
are just a change of scale of each others. In the literature there are mainly 
three typical choices, see Table [H each one with its own reason. 

Another way to obtain ([1]) is to take the random variables, Hi i ~ A/'(0, N) 
for i = 1,...,A^, and Re(if,j) ~ A^(0,iV/2), lm{Hij) ~ a/(0,A^/2) for 
l<i<j<N as independent random variables. 

For the class of real symmetric (resp. quaternionic) , one defines the 
Gaussian Orthogonal Ensemble (GOE) (resp. Gaussian Symplectic Ensemble 
(GSE)) as in Definition [1] but with [5 = 1 (resp. /9 = 4) and, of course, the 
reference measure is now the product Lebesgue measure over the independent 
entries of the matrices. 

1.2 Eigenvalues' distribution 

For the Gaussian ensembles of random matrices, one quantity of interest are 
the eigenvalues, which are independent of the choice of basis. By integrating 
out the angular variables one obtains the following result. Denote by -Pgue(A) 
the probability density of eigenvalues at A G M^. 
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Proposition 2. Let Ai, A2, • • • , Xn ^ denote the N eigenvalues of a random 
matrix H with measure (Qp. Then, the joint distribution of eigenvalues is 
given by 

PGMX) = -^\^N{\)?\{e^v(-^>?X wtthp = 2, (3) 

^ i=l ^ ' 

A7v(A) := ni<i<j<Ar(Ai ~ \) is the Vandermonde determinant, and Z^q is a 
normalization constant. 

The following identity shows why A^v is called a determinant 

A^(A) = det[Ar^,<,,.<^. (4) 

Notice that -Pgue(# e.v. G [x,x + dx]) ~ (dx)^, so that the probability of 
having double points is zero: it is a simple point process. 

For GOE (resp. GSE) the joint distributions of eigenvalues has the 
form (|3]) but with /3 = 1 (resp. /3 = 4) instead. 



1.3 Orthogonal polynomials 

Here we describe orthogonal polynomials on M, but the formulas can be 
easily adapted for polynomials on Z by replacing the Lebesgue measure by 
the counting measure and integrals by sums. 

Definition 3. Given a weight w : M 1— )■ M_|_, the orthogonal polynomials 

{gjt(a;), /c > 0} are defined by the following two conditions, 

1. qk{x) is a polynomial of degree k with qk{x) = UkX^ + . . ., Uk > 

2. they satisfy the orthonormality condition, 

{qk,qi)w-= / i^{x)qk{x)qi{x)dx = 6k,i. (5) 

Remark 4. There are other normalizations which are often used, like in 
the Askey Scheme of hypergeometric orthogonal polynomials [36]. Some- 
times, the polynomials are taken to be monic, i.e., Uk = 1 and the or- 
thonormality condition has then to be replaced by an orthogonality condi- 
tion J^u{x)qk{x)qi{x)dx = CkSk,i. Of course qk{x) = qk{x)/uk and Ck = l/ul- 
Sometimes, the polynomials are neither orthonormal (like in Definition [3]) 
nor monic, like the standard Hermite polynomials that we will encounter, 
and are given by derivatives of a given (generating) function. 
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A useful formula for sums of orthogonal polynomials is the Christoffel- 
Darboux formula: 

( UN-i qN{x)qN-i{y) - qN~i{x)qN{y) ioi x y 

qk{x)qk{y) = < ^ ~ ^ 

k=o JLJ.(^q' (^x)qN-i{x) - q'N^i{x)qN{x)), for x = y. 

K UN 

(6) 

This formula is proven by employing the following three term relation 

qn{x) = {AnX + Bn)qn~l{x) - Cnqn^2{x), (7) 

with An > 0,Bn,Cn > some constants. See Appendix [Bl for details of the 
derivation. For polynomials as in Definition [3l it holds An = Un/un-i and 

Cn = An/An-l = UnUn-2/Un-l- 

1.4 Correlation functions of GUE 

Now we restrict the point of view to the GUE ensemble and discuss the 
derivation of the correlation functions for the GUE eigenvalues' point process. 

Let the reference measure be Lebesgue. Then, the n-point correlation 
function, Pgue(^1; • • • ; ^n) is the probability density of finding an eigenvalue 
at each of the Xk, k = 1, . . . ,n. Pqve defined in ([3]) is symmetric in M^, 
which directly implies the following result. 

Lemma 5. The n-point correlation function for GUE eigenvalues is given 
by 

(n) N\ f 
PgUe(^1' • • • ' ^n) = TTT rr / PgVe{Xi, ■ ■ ■ ,XN)dXn+l ■ ■ ■ dXN (8) 

[I\ — n)\ J^N-n 

for n = 1, . . . , N and Pq\j-e{xi, ■ ■ ■ , Xn) = for n > N. 

It is important to notice that it is not fixed which eigenvalue is at which 
position. In particular, we have 

pg>,(.)= eigenvalues- dan..ty at X ^ / pgi,(x)dx = iV. (9) 

Jr 

More generally, 

[ PGUE(a^i. • • • ) Xn)dxi ...dxn= , ' . ■ . (10) 
6 



Our next goal is to do the integration in For that purpose, we 

will use Hermite orthogonal polynomials. For any family of polynomials 
{qk, k = 0, . . . , N — 1}, where qk has degree k, by multi-linearity of the de- 
terminant, it holds 

AAr(A) = det[Xl~\<ij<N = const x det[qj^i{Xi)]i<ij<N- (H) 

Therefore, setting ijj{x) := exp(— a;^/2A^), we have 

-Pgue(Ai, . . . , Xn) 

N 

= const X det[gA:-i(Ai)]i<j,fc<7vdet[gfe_i(Aj)]i<fc,j<Ar]^w(Ai) 



i=l 



12) 



const X det 



N 



^gA:-l(Ai)gfc-l(Aj 



k=l J l<iJ<N «=1 



N 



Notice that until this point, the polynomials g's do not have to be orthogonal. 
However, if we choose the polynomials orthogonal with respect to the weight 
u, then the integrations in becomes particularly simple and nice. 

Proposition 6. Let qk be orthogonal polynomials with respect to the weight 
u){x) = exp{—x'^/2N). Then, 

pSe(^i, ■■■,^n)= det x,)] , (13) 

where 

N-l 

Kf^'^ix^y) = ^^x)^:J{y)Y,(lk{x)qk{y). (14) 



A;=0 



The proof of Proposition [6] can be found in Appendix O What we need 
to do is integrate out the N — n variables one after the other and see that 
the determinant keeps the same entries but becomes smaller. For that we 
use the following two identities 



Kf^^{x,x)dx = N, 

(15) 

/ K'^''^{x,z)K^'^^{z,y)dz = K'^^^{x,y). 
Jr 
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It is at this point that the particular choice of orthogonal polynomials is 
essential. Indeed, (HM hold for the kernel defined by the r.h.s. of ([H]) precisely 
when qk{x) are taken to be orthogonal to u{x). 

The particular form of correlation functions f fT3|) appears also in other 
systems. 

Definition 7. A point process (i.e., a random point measure) is called de- 
terminantal if its n-point correlation function has the form 

p("^(a;i, ...,Xn) = det[K{xi,Xj)]i<ij<n (16) 

for some (measurable) function i^T : — )■ M, called the kernel of the deter- 
minantal point process. 

One might ask when a measure defines a determinantal point process. A 
sufficient condition is the following (see Proposition 2.2 of [1]). 

Theorem 8. Consider a probability measure on of the form 

1 ^ 
— det [0i (xj )] i<i J < AT det [V'j (xjO] i<j j< Af TT w {xi)<lxi , (17) 

1=1 



with the normalization Zj^ 7^ 0. Then (1/7) defines a determinantal point 
process with kernel 

N 

KN{x,y) = Ux)[A-%,A{y)^ (18) 

where A = [A,j]i<ij<Ar, 

Aij = = / u{z)(j)i{z)4jj{z)dz. (19) 



1.5 GUE kernel and Hermite polynomials 

Finally, let us give an explicit formula for the kernel Kn- The standard 
Hermite polynomials, {Hn,n > 0}, are defined by 



H,{y) = {-!)' ey'—e-y\ (20) 
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They satisfy 

/ e-y"Hk{y)Hi{y)dy = V^2''k\Sk,i, (21) 

JR 

with Hk{y) = 2^y^ + . . ., and also 

^ [e-y"Hr,{y)) = -e-y" Hr,^,{y) £ e'y" H^^,{y)dy = -e-^"H^{x). 

^ ~^ (22) 

By the change of variable y = x/y/^N and a simple computation, one 
shows that 

- ^vm"^ ( vif ) '''' 

are orthogonal polynomials with respect to u){x) = exp{—x'^/2N), and that 
Uk = {2TiN)-^'^k\~^''^N-^/^. Then, Christoffel-Darboux formula (ED gives 

qN{x)qN-i{y) - qN-i{x)qN{y) ^ _^,2^^2y^^ for X ^ w 

Kf^''{x^y) = { ^-y 

{qMqN-iix) - gJv-i(a:)g^(x))iVe-(^'+^')/^^, for x = y. 

(24) 

With our normalization in ([1]) the eigenvalues' density remains bounded 
and the largest eigenvalue is around the value 2N. Indeed, the eigenvalues' 
density at position fiN is given by 



-v/l-(/i/2)2, for /iG [-2,2], 

TT 



0, otherwise. 

(25) 

The large- iV asymptotic density in ( 12^ is called Wigner's semicircle law. 

1.6 Distribution of the largest eigenvalue: gap proba- 
bility 

How to compute the distribution of the largest eigenvalue, Amax? One uses 
the following simple relation, 

F'(Amax < s) = P(no eigenvalue lies in (s, oo)). (26) 
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This is a special case of gap probability, i.e., probability that in a Borel 
set B there are no eigenvalues. The gap probabilities are expressed in terms 
of n-point correlation functions as follows. 



P(no eigenvalue lies in S) = El - IsiXi)) 

^ i 

n>0 ^ ii<...<i„ k=l ' n>0 ^ n,.--,«n fc=l 

all different 

n>0 ^ 

(27) 

where 1^(2;) = 1 if x G i? and 1b (2;) = if x ^ S. The last step holds for 
simple point processes. 

For a determinantal point process, like the GUE eigenvalues, the r.h.s. of 
(I27|) is the Fredholm series expansion of the Fredholm determinant: 

P(Amax <S) = ^ tlL I det[ir^U^(Xi, X,)]i<ij<„dXi . . . dx„ 

n J{s,Oo)'^ (28) 



n=0 



= det(l — )L2((s_oo)_da:)- 

In our case the sums over n are actually finite (for n > N the terms are 
equal to zero, since the kernel Kn has rank A^), but we preferred to keep the 
formulation of the general case. 



1.7 Correlation functions of GUE minors: interlacing 
structure 

Interestingly, the determinantal structure still holds in the enlarged setting 
of eigenvalues of minors. 

Consider a N x N GUE random matrix H and denote by , . . . , Aj^ 
its eigenvalues. Then, denote by Hm the m x m minor of the matrix H 
where the last N — m rows and columns have been deleted. Denote by 
A^, . . . , the eigenvalues of Hm- In [26l[33] the correlation functions of 
{A^,l<fc<m< A^} are determined. It turns out that also in that case 
the correlation functions are determinantal (on {1, . . . , A^} x M). 



10 



4 
3 
2 
1 



n 



^2 



^2 



^^3 



»A2 



•A^ 



•A^ 



Figure 1: Illustration of the interlacing structure of the GUE minors' eigen- 
values. 



Let us order the eigenvalues for each minors as follows, 
A^ < A^ < . . . < A^. Then, the GUE minor measure can be written 
as, see e.g. 



const X ( J] det[0(Ar, Af +^)]i<,,,<™+i J det[*j^_,(Af )]i<i,,<^ (29) 

^ m=l ^ 

where A^.,.;^ = virt are virtual variables, (p{x,y) = l{x > y), 0(virt,y) = 1, 
and 

K-ki^) = (^^^)(;_,)/, e-V^^g^,,(x/y2iV). (30) 

The eigenvalues of the minors have a nice structure which appears in other 
problems too, see next sections. (!29|) is non-zero if and only if the eigenvalues 
satisfy the following interlacing structure, see Figured] for an illustration, 

Ar' < Ar < XT' < A^ < . . . < A™ < A™+}, (31) 
for m = l,...,A^— 1. Indeed, 

deti,(Ar.Ar')i..,...,...^{;: ifjr'"^'' (32) 

One might have seen that some inequalities were replaced by strict inequal- 
ities, but this is irrelevant since the events with A^ = A^"*"^ have probability 
zero. 

Let us define, 

K-ki^)--=i<P*KXlk)i^)- (33) 
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Then, a straightforward consequence of 

( 1 \n—k 

^^.(^) = jk^2^-'''''''H^-^(^/^) (34) 

for 1 < A; < n. Moreover, it is easy to verify that 



/ 



n-l 



l<k<m<n—l m=l 



(35) 

To get an exphcit expression of the kernel (see Theorem [TO] below) . we need 
to find functions k = 1, . . . ,n} orthogonal, with respect to the weight 

w = 1, to the functions j = 1, . . . , n}, and such that 

span{<l>'o^(a;), . . . , <^>^l(x)} = span{l, x,..., x^-^}. (36) 

We find 

(^l)n-j / ]\[\ ("-i)/2 



^:_^{x) = ( - ) i/„_,(x/V2iV). (37) 



2vr(n-j)! 

A measure of the form (|29l) has determinantal correlations [9J . The differ- 
ence with the above case is that the correlation functions are determinantal 
on {1, . . . , iV} X M instead of M. This means the following. The probability 
density of finding an eigenvalue of at position Xj, for i = 1, . . . , n, is given 
by 

p(")((n„ Xi),l<i<n) = det [ir^u^(n„ x,; n,-, x,)] . (38) 

The kernel is computed applying the general Theorem [10] which we 

report below, but first we give the solution. 

Proposition 9. The correlation functions of the GUE minors are determi- 
nantal with kernel given by 

Rf'^'in,, xi; n2, x,) = -(</,*(— i)) (a;,, X2)lK<n.] + ^ KUi^iWZ-ki^^) ^ 

k=l 

(39) 

with 0*'^"^ the convolution of with itself n2 — Ui times, namely, for 
n2 > ni, 

(0*("-"^))(Xi,X2) = Y l[.,-.,>0]- (40) 

177,2 — ni — i \ 
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The following result is generalization of Theorem El 

Theorem 10. Assume we have a measure on {x", 1 < i < n < N} given in 
the form, 

1 /^"' \ 

n det[Mx7,x]^%<ij<n+ij det[^^_,{xf)]i<,,,<N, (41) 

^ ^ n=l ' 

where x'^^^ are some "virtual" variables and Zfq is a normalization constant. 
If Z]\f 7^ 0, then the correlation functions are determinantal. 

To write down the kernel we need to introduce some notations. Define 

J,{ni,n2)/^ N f *---*0n2-l)(a;,l/), < na, 

\ 0, ni>n2, ^ ^ 

where (a * b) (x, y) = a{x, z)b{z, y)dz, and, for 1 < n < N , 

^L,(^):=(0^"'^^ J = l,2,...,iV. (43) 

Set (t>o{x1, x) = 1. Then the functions 

{(00 * 0^''"^)(x?,x), . . . , (0„_2 * 0("-^'"))(x:::?,x),0„_i(xr\x)} (44) 

are linearly independent and generate the n-dimensional space Vn. Define a 
set of functions {$"(x), j = 0, . . . , n — 1} spanning Vn defined by the orthog- 
onality relations 

[ <l>nx)^]ix)dx = 6,,, (45) 
Jr 

for < i, j < n — 1. 

Under Assumption (A): (f)n{Xn+i,x) = Cn^o^^^\x), for some Cn 7^ 0, 
n = l,...,A^ — 1, the kernel takes the simple form 

K{n,, xi; n2, x^) = (xi, x^) + J] KUixi)Kl-kix2). (46) 

fc=i 

Remark 11. Without Assumption (A), the correlations functions are still 
determinantal but the formula is modified as follows. Let M be the N x N 
dimensional matrix defined by [M]ij = (0j_i * * '^^_j){x\~^). Then 

K{ni,Xi]n2,X2) (47) 

N n2 

= -0("^'"^n^l, ^2) + "^nUiXl) X^[M-l],,K0,_i * 0(^-"^))(xi-\ X2). 

k=l 1=1 
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The "virtual variables" are just there to write the formula in a simpler way, 
but they do not represent real variables. This can be seen for example in 
Assumption (A), where the r.h.s. does not depend on x'^^^. Theorem [TO] is 
proven by using the framework of [TTj . 

In the case of the measure (I29p . the n-dimensional space Vn is spanned 
by {1, X, ... , x""^}, so $^ have to be polynomials of degree fc, compare with 

iZD. 

2 Totally Asymmetric Simple Exclusion Pro- 
cess (TASEP) 

2.1 Continuous time TASEP: interlacing structure 

The totally asymmetric simple exclusion process (TASEP) is one of the sim- 
plest interacting stochastic particle systems. It consists of particles on the 
lattice of integers, Z, with at most one particle at each site (exclusion prin- 
ciple). The dynamics in continuous time is as follows. Particles jump on the 
neighboring right site with rate 1 provided that the site is empty. This means 
that jumps are independent of each other and take place after an exponential 
waiting time with mean 1, which is counted from the time instant when the 
right neighbor site is empty. 

Here we consider all particles with equal rate 1. However, the framework 
which we explain below, can be generalized to particle-dependent rates and 
also particles can jump on both directions (but with different rules, it is 
not the partially asymmetric case): by jumping on their right, particles can 
be blocked, while on the left if a site is occupied, then it is pushed by the 
jumping particle. This generalization, called PushASEP, together with a 
partial extension to space-time correlations is the content of our paper [6]. 

On the macroscopic level the particle density, u{x, t), evolves deterministi- 
cally according to the Burgers equation dtu + dx{u{l—u)) = |49j. Therefore 
a natural question is to focus on fluctuation properties, which exhibit rather 
unexpected features. The asymptotic results can be found in the literature, 
see Appendix |X1 Here we focus on a method which can be used to analyze 
the joint distributions of particles' positions. This method is based on a in- 
terlacing structure first discovered by Sasamoto in |51], later extended and 
generalized in a series of papers, starting with [5]. We explain the key steps 
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following the notations of [0], where the details of the proofs can be found. 

Consider the TASEP with particles starting at time t = at positions 
Un < ■ ■ ■ < y2 < Hi- The first step is to obtain the probability that at time t 
these particles are at positions xn < ■ ■ ■ < X2 < xi, which we denote by 

. . . ,XAr;t) = P((xAr, . . . , Xi; t) | (t/at, . . . ,i/i;0)). (48) 

This function has firstly been determined using Bethe-Ansatz method [53] . A 
posteriori, the result can be checked directly by writing the master equation. 

Lemma 12. The transition probability is given by 

G{xi, . . . ,XN;t) = det{Fi_j{xN+i~i - yN+i-j,'t))i<i,j<N (49) 

with 

Ki.,t)^i^i ii^il^e'.-). (50) 

where Tq i is any simple loop oriented anticlockwise which includes w = 
and w = 1. 

The key property of Sasamoto's decomposition is the following relation 

Fn+lix,t) = Y,Fniy,t). (51) 
y>x 

Denote '■= Xk the position of TASEP particles. Using the multi- linearity 
of the determinant and (l5Ti) one obtains 

G{xi, . . . , xtv; t) = ^ det(F_j+i(xf - yN-j+i, t))i<ij<N, (52) 
v 

where 

V' = {x^, 2<i<n< N\x1 > x^Zi}. (53) 
Then, using the antisymmetry of the determinant and Lemma [13] below we 



can rewrite (I52|) as 

G{xi, . . . , xat; t) = ^ det(F_j+i(xf - yN-j+i, t))i<i,j<N, (54) 



where 



V = {x1, 2<i<n< N\x'l > <+\ > xti}. (55) 
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Lemma 13. Let f an antisymmetric function of {x^ , . . . ,x^}. Then, when- 
ever f has enough decay to make the sums finite, 

5^/(xf,...,4) = 5^/(xf,...,4), (56) 
V v 

with the positions x\> x\> . . . > x^ being fixed. 

Now, notice that, for n = —k < 0, (|5CT]) has only a pole at w = 0, which 
implies that 

Fn+iix,t) = -Y,Fn{x,t). (57) 

y<x 

Define then 

1 f dw (1 — w^'^ 



2m /r^ w w^'-y^-"-^ (58) 
0(x, y) := l{x > y), 0(virt, y) = 1. 
In particular, for /c = 1, . . . , iV, (x) = { — l)^F^k{x — yN-k, t)- Then, 

where the interlacing condition V, is given by 

< < <x-,<...<x:< x^Xl, (60) 

for 77, = 1, . . . , — 1, see Figure [2] for a graphical representation. 

Therefore, we can replace the sum over I? by a product of determinants 
of increasing sizes. Namely, 

G(xi,...,x^;t) = Qi{xkA<k<n<N}) (61) 

2<k<n<N 

where the measure Q is 

N-l 

Qi{xl l<k<n<N}) = (^l[ det(0(x^x^"+i))i<,,,<„+i' 

n=l 

N 

where we set x^_^_i = virt (we call them virtual variables). As for the GUE 
minor case, if we compute \E'^_^(x) := (0 * ^n+i-fe)!^)? S^t the fl58|) but 
with replaced by n. 



(62) 



X det(^j^_,(xf ))i<i,j<jv, 
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Figure 2: Graphical representation of the domain of integration T) for = 4. 
One has to "integrate" out the variables xl, i>2 (i.e., the black dots). The 
positions of xj, = 1, . . . , are given (i.e., the white dots). 



Remark 14. The measure f l62|) is not necessarily a probability measure, since 
positivity is not ensured, but still has the same structure as the GUE minors 
measure, compare with f l2^ . In particular, the distribution of TASEP parti- 
cles' position can be expressed in the same form as the gap probability, but 
in general for a signed determinantal point process. The precise statement 
is in Theorem [15] below. 

Applying the discrete version of Theorem [TU| see Lemma 3.4 of [S], we 
get the following result. 

Theorem 15. Let us start at time t = with N particles at positions 
Un < ■ ■ ■ < y2 < yi- Let cr(l) < cr(2) < . . . < a{m) he the indices of m 
out of the N particles. The joint distribution of their positions Xfj{k){t) is 
given by 

m 

P( Pi {Xa(k){t) > Sk}^ = det(l - XsKtXs) email),.. .,airn)}xZ) (63) 
k=l 

where Xs{o'{k),x) = l(x < Sk). Kt is the extended kernel with entries 

"2 

Kt{n,, xi; n2, x,) = -0("-"^)(a;i, x,) + ^ KUMKUM (64) 

k=l 

where 

0(---)(xi,X2)= f'^^^'^^^lY (65) 
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and the functions ■i = 0,...,?T, — 1, /orm a family of polynomials of 

degree < n satisfying 

J2mx)<^li^)=kr (67) 

The path Tq in the definition of'^f is any simple loop, anticlockwise oriented, 
which includes the pole at w = but not the one at w = 1. 

We omit to write explicitly the dependence on the set {yi} is hidden in 
the definition of the $"'s and the ^'"'s. 

2.2 Correlation functions for step initial conditions: 
Charlier polynomials 

Now we consider the particular case of step initial conditions, yk = —k for 
k > 1. In that case, the measure ( l62l) is positive, i.e., it is a probability 
measure. The correlation functions of subsets of {x^{t), 1 < k < n,n > 1} 
are determinantal with kernel Kj^^^^, which is computed by Theorem [T5l 

The correlation kernel /^tasep gjyg^ terms of Charlier polynomials, 
which we now introduce. The Charlier polynomial of degree n is denoted by 
Cn{x, t) and defined as follows. Consider the weight a; on Z+ = {0, 1, . . .} 

uj{x) = e'Hyxl. (68) 

Then, Cn{x,t) is defined via the orthogonal condition 

V — V Tl\ 

2_^Uj{x)Cn{x,t)Cm{x,t) = —6n,m (69) 

a;>0 

or, equivalently, C„(x, t) = (— l/t)"x" H — ■ . They can be expressed in terms 
of hypergeometric functions 

Cn{x, t) = 2Fo(-n, -x; ; -1/t) = C,{n, t). (70) 

From the generating function of the Charlier polynomials 

^ Cnix^^^^^n ^ ^^t^^ _ (7;l) 
n>0 ^' 
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one gets the integral representation 



-,Cn{x,t) = ^f dw '^'^\j^\ (72) 
n! zvri 

Equations ( ITOl) and ( 172|) then give 

^^(-n + a;) = ^^Ck{x, t) and ^^"(-n + a;) = jCj{x, t). (73) 
In particular, the kernel of the joint distributions of {x^ , . . . , x^} is 

N-l 



^TASEP(^^ -AT + x; iV, -iV + y) = ^ ^^{-N + x)$^(-iV + y) 

N-l 

= uj{x)^qk{x)qk{y) 



(74) 



k=0 



with gfc(a;) = (— l)^^Cfc(x, t) = {t^B?)~^/'^x^ + . . ., orthogonal polynomials 
with respect to uj{x). Therefore, using the Christoffel-Darboux formula ([6]) 
we get, for x ^ y, 



^TASEP^^^ -A^ + x; AT, -A^ + y) 

^ e-H- CN-iix,t)CN{y,t)-CN{x,t)CN-iiy,t) 
x\ (A^-1)! x-y 



■ (75) 



Remark 16. The expression of the kernel in (17^ is not symmetric in x 
and y. Is this wrong? No! For a determinantal point process the ker- 
nel is not uniquely defined, but only up to conjugation. For any function 
/ > 0, the kernel K{x,y) = f (x) K {x , y) f {y)~^ describes the same determi- 
nantal point process of the kernel K{x,y). Indeed, the / cancels out in the 
determinants defining the correlation functions. In particular, by choosing 
f{x) = 1/ >/uj{x), we get a symmetric version of (j74l) . 

Double integral representation of /^tasep 

Another typical way of representing the kernel /('TASep jg ^^Q^i^ig integral. 
This representation is well adapted to large-t asymptotic analysis (both in 
the bulk and the edge). 
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Let us start with fl^ with = —k. We have 

and remark that ^E'^(x) = for x < —n, k > 0. The orthogonal functions 
k = 0, ... ,n — 1, have to be polynomials of degree k and to satisfy 
the orthogonal relation ^^.>_„ ^^(a;)^^ (a;) = Sij. They are given by 

_1 r p-t{z-l) x+n 

= — d> dz- (77) 

Indeed, for any choice of the paths rojFi such that \z\ < \w\, 

-1 f , f , e*('^-i) (1 - 1 
dz (h aw- 



(27ri)2 Jr^ J^^^ e*(^-i) (1 - w-z 



27ri 



(78) 

since the only pole in the last w-integral is the simple pole at w = z. 

After the orthogonalization, we can determine the kernel. Let us compute 
the last term in 1^^. For k > n2, we have $^2_^(x) = 0. So, by choosing z 
close enough to 1, we have |1 — -z] < |1— and we can extend the sum over 
k to +00, i.e., 

E = E p-p d„, f^^ ^ _ 

dw (p dz- -— — > 



(27ri)2 Jp^ /p^ ete^a;i+ni+l ^ (1 - ;z)"2-A:+l 

dw (p dz- — — > 

/„ ptz,,,,xi+ni+l 



(27ri)2 /p^^ gfey^aii+ni+l ZL-' (1 



^ dwi dz- ^ 



(27ri)2 Jp^ Jp^ e*^(l - z)"2y;xi+ni+i ;2 _ • 

(79) 
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The integral over Fi contains only the pole z = 1, while the integral over Tq 
only the pole w = (i.e., -u; = 2; is not inside the integration paths). This is 
the kernel for rii > n2- For rii < n2, there is the extra term — ^("^'"■^^(xi, X2)- 
It is not difficult to check that 



az- 



(27ri)2 Jp^ Jp^ e*^(l - ^)"2w^i+'^i+i 2: - w 
^ dw- — ' 



27ri Jp^ -u;xi+m-(3;2+n2)+i ^2 - ni - 1 

Therefore, the double integral representation of Kj^^^^ (for step initial con- 
dition) is the following: 



(27ri)2 /p^ "^^6*^(1 - ^)"2^yxi+ni+l _2 _ 

(27ri)2 fp^ ^ "~ e*^(l - 2;)'^2^i;Xi+ni+l 2; _ 



1 / / et«'(l _-u;)ni^X2+n2 1 

e*"'(l - w)"lZ^2+n2 1 
d^^TT^ .,^.,^1 > if ril < n2. 



(81) 



Remark 17. Notice that in the GUE case, the most natural object are the 
eigenvalues for the N x N matrix, , . . . , X^. They are directly associated 
with a determinantal point process. The corresponding quantities in terms 
of TASEP are not the positions of the particles, xl,...,Xi. The measure on 
these particle positions is not determinantal, but with the extension to the 
larger picture, namely to {x"^, 1 < k < n < N}, we recover a determinantal 
structure. This is used to determine the joint distributions of our particles, 
since in terms of 1 < k < n < N} we need only to compute a gap 
probability. 

2.3 Discrete time TASEP 

There are several discrete time dynamics of TASEP from which the contin- 
uous time limit can be obtained. The most common dynamics are: 

• Pcirallel update: at time t e Z one first selects which particles can 
jump (their right neighboring site is empty). Then, the configuration 
of particles at time t + 1 is obtained by moving independently with 
probability p e (0, 1) the selected particles. 
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• Sequential update: one updates the particles from right to left. The 
configuration at time t + 1 is obtained by moving with probability 
p G (0, 1) the particles whose right rite is empty. This procedure is 
from right to left, which implies that also a block of m particles can 
move in one time-step with probability p"^. 

Other dynamical rules have also been introduced, see [54j for a review. 
Sequential update 

For the TASEP with sequential update, there is an analogue of Lemma [121 
with the only difference lying in the functions (see |1H] and Lemma 3.1 
of [H]). The functions F„ satisfy again the recursion relation (I^Tj) and The- 
orem [15] still holds with the only difference being in the \l/"(x)'s which are 
now given by 



For step initial conditions, the kernel is then given by ( 18T1) with e^^/e^^ re- 
placed by {1 — p + pw f/ {1 — p + pz f. 

Parallel update 

For the TASEP with parallel update, the same formalism used above can still 
be applied. However, the interlacing condition and the transition functions 
0„ are different. The details can be found in Section 3 of [TO] (in that paper 
we also consider the case of different times, which can be used for example 
to study the tagged particle problem). The analogue of (j62l) is the following 
(see Proposition 7 of [TO]). 

Lemma 18. The transition probability G{x; t) can be written as a sum over 




(82) 



V" = {x1, 2<i<n< N\x1 > J 



(83) 



as follows: 




(84) 
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where 

N~l 

Q(K, l<k<n<N}) = det(0»(xr_i,x^"+i))i<,,,<„+i) 

n=l 

X det(F_j+i(xf - VN-j+i, t+1- j))i<i,j<N- 

(85) 

where we set Xq = — oo (we call them virtual variables). The function (p^ is 
defined by 

{1, y>x, 
p, y = x-l (86) 
0, y<x-2, 

and F_n is given by 

If 



The product of determinants in ( ]85l) also implies a weighted interlacing 
condition, different from the one of TASEP. More precisely, 

^n+i Kx'l-IK (88) 

The difference between the continuous-time TASEP is that it can happens 
that x^^^ = — 1. However, for each occurrence of such a configuration, the 
weight is multiplied by p. The continuous-time limit is obtained by replacing 
t by t/p and letting p — )■ 0. 

Also in this case Theorem [15] holds with the following new functions. 

If 1 f w x^i-'^a 



27ri Jro_i (1 + w)^-y"->+i \(1 + u;)(l -Fpw) 

(89) 

and $^ are polynomials of degree k given by the orthogonality condition (1671) : 
^xez (^) = ^id- particular, for step initial conditions, yk = —k, 

k > 1, we obtain the following result. 

Proposition 19. The correlation kernel for discrete-time TASEP with par- 
allel update is given by 

KTAS^P(ni,Xi;n2,X2) = -^^^-'^^Ha^i, 2:2) + i^i^^'^^(ni, Xi; ns, Xs), (90) 
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with (f)^' 



ni,n2) 



given by 



[n2>ni] 



i 



dw 



(1 +pw)"2-"it(;"i-"2 



(91) 



27ri 




and 




(ni,Xi;n2,X2) 



(92) 



Remark 20. The discrete-time parallel update TASEP with step initial 
condition is equivalent to the shuffling algorithm on the Aztec dynamics 
as shown in |?T]. This particle dynamics also fits in the framework devel- 
oped in [S]. See [21] for an animation, where the particles have coordinates 
(zf := + n,n). 

3 2 + 1 dynamics: connection to random 
tilings and random matrices 

In recent years there has been a lot of progress in understanding large time 
fluctuations of driven interacting particle systems on the one-dimensional 
lattice. Evolution of such systems is commonly interpreted as random growth 
of a one-dimensional interface, and if one views the time as an extra variable, 
the evolution produces a random surface (see e.g. Figure 4.5 in [44j for a nice 
illustration). In a different direction, substantial progress have also been 
achieved in studying the asymptotics of random surfaces arising from dimers 
on planar bipartite graphs. 

Although random surfaces of these two kinds were shown to share certain 
asymptotic properties (also common to random matrix models), no direct 
connection between them was known. We present a class of two-dimensional 
random growth models (that is, the principal object is a randomly growing 
surface, embedded in the four- dimensional space-time). 

In two different projections these models yield random surfaces of the two 
kinds mentioned above (one reduces the spatial dimension by one, the second 
projection is to fixced time). 

Let us now we explain the 2 + 1-dimensional dynamics. Consider the set 
of variables {x^(t), 1 < < n, > 1} and let us see what is their evolution 
inherited from the TASEP dynamics. 
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3.1 2 + 1 dynamics for continuous time TASEP 



Consider continuous time TASEP with step-initial condition, yu = —k for 
k>l. Then, a further property of the measure fl62|) with step initial condi- 
tions is that 

xl{Q) = -n + k-l. (93) 

Let us verify it. The first N — 1 determinants in fl62|) imply the interlacing 
condition ( 160|1 . In particular, x]^(0) = — and (0) > — + k — 1 for 
k > 2. At time t = we have 

FJx,0) = — i dw ^"" ~}}^ . (94) 

By residue's theorem, we have F_fc(x, 0) = for x > 1 (no pole at oo), 
F_fc(x, 0) = for X < -k (no pole at 0), and F_fc(0, 0) = 1. 

The last determinant in fl62p is then given by the determinant of 

Fo(0,0) F-i(-l,0) ■■■ F_^+i(-Ar + l,0) 

Fo(x^(0) + Ar,0) F_i(x^(0) + Ar-1) ... F„;v+i(x^(0) + 1, 0) 

Fo(xj^(0) + Ar,0) F_i(xj^(0) + Ar-1,0) ... F_;v+i(xj^(0) + 1, 0), 

(95) 

Let us see when det( l95|) ^ 0: 

1. since (0) > -A^ + A; - 1, the first column of (ESI) is [1, 0, ... , 0]*. 

2. Then, if x^(0) > — A^ + 1, the second column is [*,0,...,0]* and 
det(|95|) = 0. Thus we have x^(0) = -A^ + 1. 

3. Repeating the argument for the other columns, we obtain that 
detdlSD ^ if and only if (0) = -A^ + A; - 1 for = 3, . . . , A^. 

This initial condition is illustrated in Figure [3] (top, left). 

Now we explain the dynamics on the variables {x^(t), 1 < /c < n,n > 1} 
which is inherited by the dynamics on the TASEP particles {x"(t),n > 1}. 
Each of the particles x^ has an independent exponential clock of rate one, 
and when the x™-clock rings the particle attempts to jump to the right by 
one. If at that moment x™ = x^~^ — 1 then the jump is blocked. If that is 
not the case, we find the largest c > 1 such that x^ = x^^^ = ■ ■ ■ = x^'i'i^, 
and all c particles in this string jump to the right by one. 
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Figure 3: (top, left) Illustration of the initial conditions for the particles 
system, (bottom, left) A configuration obtained from the initial conditions, 
(right) The corresponding lozenge tiling configurations. In the height func- 
tion picture, the white circle has coordinates {x,n,h) = (—1/2,0,0). For a 
Java animation of the model see pU] . 
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Informally speaking, the particles with smaller upper indices are heavier 
than those with larger upper indices, so that the heavier particles block and 
push the lighter ones in order for the interlacing conditions to be preserved. 

We illustrate the dynamics using Figure |3l which shows a possible con- 
figuration of particles obtained from our initial condition. If in this state of 
the system the x^-clock rings, then particle xf does not move, because it is 
blocked by particle xf. If it is the a;|-clock that rings, then particle moves 
to the right by one unit, but to keep the interlacing property satisfied, also 
particles x| and xf move by one unit at the same time. This aspect of the 
dynamics is called "pushing". 

3.2 Interface growth interpretation 

Figure |3] (right) has a clear three-dimensional connotation. Given the random 
configuration {x'^{t)} at time moment t, define the random height function 

/i : (Z + i) X Z>o X M>o ^ Z>o, 

h{x,n,t) = 4^{ke{l,...,n}\xUt)>x}. ^ ^ 

In terms of the tiling on Figure |3l the height function is defined at the vertices 
of rhombi, and it counts the number of particles to the right from a given 
vertex. (This definition differs by a simple linear function of {x, n) from the 
standard definition of the height function for lozenge tilings, see e.g. [Mt[35].) 
The initial condition corresponds to starting with perfectly fiat facets. 

In terms of the step surface of Figure [31 the evolution consists of re- 
moving all columns of (x, n, /i)-dimensions (1,*,1) that could be removed, 
independently with exponential waiting times of mean one. For example, if 
X2 jumps to its right, then three consecutive cubes (associated to x1,x\,xX) 
are removed. Clearly, in this dynamics the directions x and n do not play 
symmetric roles. Indeed, this model belongs to the 2 + 1 anisotropic KPZ 
class of stochastic growth models, see [5l[^. 

3.3 Random tilings interpretation 

A further interpretation of our particles' system is a random tiling model. To 
see that one surrounds each particle location by a rhombus of one type (the 
light-gray in Figure [3]) and draws edges through locations where there are no 
particles. In this way we have a random tiling with three types of tiles that 
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Figure 4: Illustration of the dynamics on tiles for a column of height m = 4. 



we call white, light-gray, and dark-gray. Our initial condition corresponds to 
a perfectly regular tiling. 

The dynamics on random tilings is the following. Consider all sub- 
configuration of the random tiling which looks like a visible column, i.e., 
for some m > 1, there are m light-gray tiles on the left of m white tiles (and 
then automatically closed by a dark-gray tile). The dynamics is an exchange 
of light-gray and white tiles within the column. More precisely, for a column 
of height m, for all k = 1, . . . ,m, independently and with rate 1, there is 
an exchange between the top k light-gray tiles with the top white tiles as 
illustrated in Figure H] for the case m = 4. 

Remark 21. We can also derive a determinantal formula not only for the 
correlation of light-gray tiles, but also for the three types of tiles. This is 
explicitly stated in Theorem 5.2 of [5]. 

3.4 Diffusion scaling and relation with GUE minors 

There is an interesting partial link with GUE minors. In the diffusion scaling 
limit 



the measure on {C,]^, 1 < k < n < N} is exactly given by f l29|) . 

Remark 22. It is important to stress, that this correspondence is a fixed- 
time result. From this, a dynamical equivalence does not follow. Indeed, if 




(97) 
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Figure 5: A random tiling of the Aztec diamond of size n = 10. 



we let the GUE matrices evolves according to the so-called Dyson's Brownian 
Motion dynamics, then the evolution of the minors is not the same as the 
(properly rescaled) evolution from our 2 + 1 dynamics for TASEP. 

3.5 Shuffling algorithm and discrete time TASEP 

An Aztec diamond is a shape like the outer border of Figure |5l The shuffling 
algorithm provides a way of generating after n steps uniform sample of an 
Aztec diamond of size n tiled with dominos. 

We now discuss (not prove) the connection between discrete time TASEP 
with parallel update and step initial condition. Moreover, we take the pa- 
rameter p = 1/2 to get uniform distribution. It is helpful to do a linear 
change of variable. Instead of we use 

z^ = xl + n, (98) 

so that the interlacing condition becomes 

4"-' < 4 < 4X1 (99) 

The step initial condition for TASEP particles is Zi{Q) = 0, n > 1. An 
analysis similar to the one of Section |3TT] leads to 2:^(0) = k — l,l<k<n. 
Then, the dynamics on {z]^,l < k < n,n > 1} inherited by discrete time 
parallel update TASEP is the following. First of all, during the time-step 
from n — 1 to n, all particles with upper-index greater or equal to n + 1 are 
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t=2 t=l t=0 t=l t=2 




Figure 6: Two examples of configurations at time t = 2 obtained by the 
particle dynamics and its associated Aztec diamonds (rotated by 45 degrees). 



frozen. Then, from level n down to level 1, particles independently jump 
on their right site with probability 1/2, provided the interlacing condition 
fl99|) with the lower levels are satisfied. If the interlacing condition would be 
violated for particles in upper levels, then these particles are also pushed by 
one position to restore (1^ . 

Finally, let us explain how to associate a tiling configuration to a particle 
configuration. Up to time t = n only particles with upper-index at most n 
could have possibly moved. These are also the only particles which are taken 
into account to determine the random tiling. The tiling follows these rules, 
see Figure E] for an illustration: 

1. light-gray tiles: placed on each particle which moved in the last time- 
step, 

2. middle-gray tiles: placed on each particle which did not move in the 
last time-step, 

3. dark-gray tiles and white tiles: in the remaining position, depending 
on the tile orientation. 

The proof of the equivalence of the dynamics can be found in [H] , where 
particle positions are slightly shifted with respect to Figure [61 In [21] you 
can find a Java animation of the dynamics. 
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A Further references 



In this section we give further references, in particular, of papers based on 
the approach described in this lecture notes. 

• Interlacing structure and random matrices: In [33] it is studied the 
GUE minor process which arises also in the Aztec diamond at the 
turning points. Turning points and GUE minor process occurs also 
for some class of Young diagrams [12]. The antisymmetric version of 
the GUE minors is studied in [25] . In [26] the correlation functions for 
several random matrix ensembles are obtained, using two methods: the 
interlacing structure from |9j and the approach of PU]. When taking 
the limit into the bulk of the GUE minors one obtains the bead process, 
see [18]. Further works on interlacing structures are [T9l[32ll39| [56] . 

• GUE minors and TASEP: Both the GUE minor process and the anti- 
symmetric one occurs in the diffusion scaling limit of TASEP [TT|[T2]. 

• 2+1 dynamics: The Markov process on interlacing structure introduced 
in |5j is not restricted to continuous time TASEP, but it is much more 
general. For example, it holds for PushASEP dynamics [H] and can 
be used for growth with a wall too [15]. In a discrete setting, a similar 
approach leads to a shuffling algorithm for boxed plane partitions [T3] . 
As already mentioned, the connection between shuffling algorithm and 
interlacing particle dynamics is proven in [?T] (the connection with 
discrete time TASEP is however not mentioned). 

• 2 + 1 anisotropic growth: In the large time limit in the 2 + 1 growth 
model the Gaussian Free Field arises, see [5J or for a more physical 
description of the result [7]. In particular, height fluctuations live on a 



vhit scale (in the bulk) and our model belongs to the anisotropic KPZ 

class, like the model studied in [45] . 

• Interlacing and asymptotics of TASEP: Large time asymptotics of 
TASEP particles' positions with a few but important types of ini- 
tial condition have been worked out using the framework initiated 
with [9]. Periodic initial conditions are studied in [9] and for discrete 
time TASEP (sequential update [8], parallel update [10] )• The limit 
process of the rescaled particles' positions is the Airyi process. For step 
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initial condition it was the Airy2 process [21] • The transition process 
between these two has been discovered in [TT] , see also the review [23] . 
Finally, the above technique can be used also for non-uniform jump 
rates where shock can occur [T2| . 

• Line ensembles method and corner growth models: TASEP can be also 
interpreted as a growth model, if the occupation variables are taken to 
be the discrete gradient of an interface. TASEP belongs to the so-called 
Kardar-Parisi- Zhang (KPZ) universality class of growth models. It is 
in this context that the first connections between random matrices and 
stochastic growth models have been obtained [27j. The model studied is 
analogue to step initial conditions for TASEP. This initial condition can 
be studied using non-intersection line ensembles methods [28l[29]. The 
Airy2 process was discovered in [17] using this method, see also [3T | [52 | 
[55] for reviews on this technique. The non-intersecting line description 
is used also to prove the occurrence of the Airy2 process at the edge of 
the frozen region in the Aztec diamond |30] . 

• Stationary TASEP and directed percolation: Directed percolation for 
exponential/geometric random variables is tightly related with TASEP. 
In particular, the two-point function of stationary TASEP can be re- 
lated with a directed percolation model ^6]. The large time behavior 
of the two-point function conjectured in [46] based on universality was 
proven in [24j. Some other universality-based of [IB] conjectures have 
been verified in [T]. The large time limit process of particles' positions 
in stationary TASEP, the corresponding point-to-point directed perco- 
lation (with sources), and also for a related queuing system, has been 
unraveled in [3]. The different models shares the same asymptotics due 
to the slow-decorrelation phenomena [22] . 

• Directed percolation and random matrices: Directed percolation, the 
Schur process and random matrices have also nice connections; from 
sample covariance matrices [2], to small rank perturbation of Hermitian 
random matrices [M], and to the generalization [IB] . 
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B ChristofFel-Darboux formula 

Here we prove Christoffel-Darboux formula ([6]). First of all, we prove the 
three term relation ([7]). From = x" + ■ ■ ■ it follows that 

Qnjx) _ xqn-i{x) ^^^^^ 

is a polynomials of degree n — 1. Thus, 

qn{x) xqn-i{x) / qn \ 

= ^ } akqk{x), ak = { ,qk) , (101) 

Un Un-l ^ \Un Mn-1 

where X is the multiplication operator by x, and (/, g)ui = co{x)f{x)g{x)dx 
is the scalar product. 

Let us show that afc = for A; = 0, . . . , n — 3. Using {Xf, g)^ = (/, Xg)i^ 
we get 

«fc = — {qn,qk)uj {qn-i,Xqk)uj = (102) 

for /c + 1 < n — 1, since Xqk is a polynomial of degree A; + 1 and can be 
written as linear combination of go, • • • , ^fc+i- 
Consider next k = n — 2. We have 

an-2 = —{qn-l,Xqn^2)uj = -^^^^, (103) 

Un—l ^n—l 

because we can write 

xqn-2{x) = Un-2x"'~^ + a polynomial of degree n — 2 

= " qn~i{x) + a polynomial of degree n — 2. (^0^) 

Un-l 

Therefore, setting = a^-i^in, = Mn/^n-i, and C„ = UnUn-2/ul-i, we 
obtain the three term relation ([71). We rewrite it here for convenience, 

qn{x) = (AnX + Bn)qn-l{x) - Cnqn-2{x). (105) 

From (11051) it follows 

qn+l{.x)qn{y) - qn{.x)qn+l{y) 

= An+iqn{x)qn{y){x - y) + Cn+1 {qn{x) qn^l{y) " qn-l{x) qn{y)) ■ (106) 



33 



We now consider the 7^ y. The case x = y is obtained by taking the 

y ^ X hmit. Dividing ( I106p by (x — y)An+i we get, for /c > 1, 

qkix)qkiy) = Sk+iix, y) - Sk{x, y), (107) 

where we defined 

S,{x, y) = - q>^-Ax)%{y) _ (108) 

Uk X - y 

Therefore (for x ^ y) 

N-l 

^ qk{x)qk{y) = Sn{x, y) - Si{x, y) + qo{x)qo{y) = SNix, y). (109) 

fc=0 

The last step uses qo{x) = uq and qi{x) = uix + c (for some constant c), 
from which it follows qo{x)qo{y) = Si{x,y). This ends the derivation of the 
Christoffel-Darboux formula. 



C Proof of Proposition [6 



/ 



Here we present the details of the proof of Proposition [6] since it shows how 
the choice of the orthogonal polynomial is convenient. The basic ingredients 
of the proof of Theorem [8] are the same, with the only important difference 
that the functions in the determinants in f|T7|) are not yet biorthogonal. 
First of all, let us verify the two relations f|T5|) . We have 

N-l 

i^^u^(x,x)dx = ^(gfc,gfc). = N, (110) 

k=0 

and 

N-l 

K^''''{x,z)K^''^{z,y)dz= ^ ^u{x)u{y)qkix)qi{y){qk,qi)^ 

k,l=0 

= Kf'\x,y). 

By Lemma [5], Equation (fT2|) . and the definition of K^^^, we have 

PgUe(-^1' • • • ' ^n) 



;iii) 



^N ,J^'., [ det[K^^^{xi,Xj)]^^^ dXn+i...dxN. (112) 
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We need to integrate N — n times, each step is similar. Assume therefore 
that we aheady reduced the size of the determinant to m x m, i.e., integrated 
out Xm+i, ■ ■ ■ ,xn- Then, we need to compute 



/ det 



l<i,j<m 



dXr 



(113) 



m— 1 



In what follows we write only K instead of Kj^ . We develop the determi- 
nant along the last column and get 

det [K{xi, Xj)]^^^^.^^ = K{x^, Xm) det [K{xi, Xj)]^^..^^_^ 

[K{Xm, ^j)]l<j<m-l 
K{Xm, Xm) det [K(^Xi, 3^j)]i<j j<„_i 



J2{-ir-'K{xk,Xm) det 



k=l 



m—1 



E(-i) 



m—k 



det 



k=l 



i<j<m-l 



4) 



Finally, by using the two relations ( TT^ . Equation f lll3p becomes 

^-Ij-'-^det 

k=l 



[K{xk,Xj)] 



Ndet [K{x.,x,)],^^^^^^_, + J](-l)-^det 
= (AT - (rri - 1)) det [K{x,, . 

This result, applied for m = A^, — 1, . . . , + 1, leads to 

c^vA^ldet [K^''''{x„x,)] ^ 



l<j<m— 1 



PgUeI^I? • • • , Xn) 



(115) 



(116) 



Now we need to determine cat. Since cat depends only of A^, we can 
compute it for the n = 1 case. From the above computations, we have 
Pgue(^) — cnN\K^^^{x,x) and PgueI^)*^^ = A^ we have cat = 1/A^!. 
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